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A quantum field theory warm inflation model is presented 
that solves the horizon/flatness problems. The model obtains, 
from the elementary dynamics of particle physics, cosmolog- 
ical scale factor trajectories that begin in a radiation domi- 
nated regime, enter an inflationary regime and then smoothly 
exit back into a radiation dominated regime, with nonnegli- 
gible radiation throughout the evolution. 
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The resolution of the horizon problem, which underlies 
inflationary cosmology jlj, is that at a very early time, 
the equation of state that dictates the expansion rate of 
the Universe was dominated by a vacuum energy density 
p Vl so that a small causally connected patch could grow 
to a size that encompasses the comoving volume which 
becomes the observed universe today. 

In the standard (isentropic) inflationary scenarios, the 
radiation energy density p r scales with the inverse fourth 
power of the scale factor, becoming quickly negligible. In 
such short time reheating period terminates the 

inflationary period and initiates the radiation dominated 
epoch. On the other hand, the only condition required by 
General Relativity for inflation is that p r < p v . Inflation 
in the presence of nonnegligible radiation is characterized 
by non-isentropic expansion |^,|3| and thermal seeds of 
density perturbations This can be realized in warm 
inflation scenarios ^ where there is no reheating. 

The basic idea of our implementation of warm inflation 
is quite simple; a scalar field, which we call the inflaton, 
is coupled to several other fields. As the inflaton relaxes 
toward its minimum energy configuration, it will decay 
into lighter fields, generating an effective viscosity. That 
this indeed happens has been demonstrated in detail in 
Refs. [|| ||. If this viscosity is large enough, the inflaton 
will reach a slow-roll regime, where its dynamics becomes 
overdamped. This overdamped regime has been analyzed 
in Ref. ||. As one expects, overdamping is most success- 
ful for the case where the inflaton is coupled to a large 
number of fields which are thermally excited, i.e., which 
have small masses compared to the ambient temperature 
of radiation. This result has important consequences for 
cosmological applications. 



In order to satisfy one of the requirements of a suc- 
cessful inflation (60 or so e- folds), overdamping must be 
very efficient. Our goal in this Letter is to show that 
it is possible to build a toy model, motivated from high 
energy particle physics, that can provide enough over- 
damping as to produce a viable inflationary expansion, 
which smoothly exits into a radiation-dominated regime, 
with p r slowly and monotonically decreasing throughout 
the whole process. In contrast to most models in the lit- 
erature, warm inflation provides both a natural context 
for the slow-roll regime and for its graceful termination 
into a radiation-dominated era. 

The particle physics model considered below is inspired 
by string theories exhibiting N = 1 global supersymme- 
try, with the inflaton coupled to massive modes of the 
string, as recently discussed in Ref. pS. We refer the 
reader to this reference for details. 

We thus consider the following Lagrangian of a scalar 
field (j> interacting with Nm X N x scalar fields Xjk and 
Nm x N^p fermion fields ipjk, 

1 m 2 A 

^[<P,Xjk,tpjk,ipjk] = ^(d^) 2 - —cj} 2 - 

N M N x f 2 >| 

■J=l k=l I J 

N M -W* 

+ {^Jk ^j* ~ h jk{4'- M^tpjklpjk} , (1) 

j=l k=l 

where all coupling constants are positive: A, fjk,g 2 kjhjk 
> 0. For simplicity, we consider in the following /,•/. = /, 
9jk = = g. Also, we will set Nj, — N x /A, which 
along with our choice of coupling implies a cancellation 
of radiatively generated vacuum energy corrections in the 
effective potential |ll[] . We call this kind of model a dis- 
tributed mass model (DMM), where the interaction be- 
tween <p with the Xjk and ipjk fields establishes a mass 
scale distribution for the Xjk and ipjk fields, which is de- 
termined by the mass parameters {M}. Thus the Xjk and 
ipjk effective field-dependent masses, m Xjk (<j>, T, {M}) 
and m^ jk (0, T, {M}), respectively, can be constrained 
even when ((f) = ip is large. A specific choice of {M} will 
be given shortly. The 4>x, ^ interactions can be made 
reflection symmetric, <f> — > —(f), but for our purposes we 
will consider all Mj > and (j> > 0. 

The 1-loop effective equation of motion for the scalar 
field (f> is obtained by setting (f> = ip + r] in Eq. (|) and 
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imposing (rj) — 0. Then from Weinberg's tadpole method 
the 1-loop evolution equation for ip is 



ip + 3H<p + m 2 ip + ^tp 3 + ^ip(r, 2 ) 

N M N x N m 

+9 2 E I> - M *K4-> (V'tf^) = ■ (2) 

» j i 3 

In the above, the term 3Htp describes the energy red-shift 
of ip due to the expansion of the Universe. In the warm- 
inflation regime of interest here, the characteristic time 
scales (given by the inverse of the decay width) for the 
fields in Eq. (|l|) are faster than the expansion time scale, 
X/H. In this case, the calculation of the (renormalized) 
thermal averages in Eq. (|J) can be approximated just as 
in the Minkowski spacetime case. 

A systematic evaluation of the averages in the adia- 
batic, strong dissipative regime was presented in || and 
re-derived in (m) with extension to fermions. The es- 
sential mechanism for dissipation obtained from this ap- 
proach can be explained through an intuitive kinetic the- 
ory argument first given in Q and reexamined in p4[ . 
(We note that the objections to warm inflation raised in 
Ref. Q are avoided quite naturally by coupling the infla- 
ton to a tower of mass modes as in the present implemen- 
tation of the DM model.) For the %-field averages one 
writes (x?.) = f(d 3 fc/(27r) 3 )n Xi3 (k)/w Xij .(k), where the 
number densities n(k) in the strong dissipative regime, 
in near equilibrium, can be written in the relaxation time 
approximation to the kinetic equation as plfL5| 



n(k) 



■(k) 



• eq 



■(k). 
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Here n eq (k) = l/(e ,3 ^ k ) — 1) is the equilibrium number 
density for x particles, w(k) = \/k 2 + m 2 (ip, T, {M}), 
to 2 {ip, T, {MY) is the effective, field dependent mass for 
the %- fields §), and r = where T is the decay width 
for the ^-particles. From Eq. (^) the second term on 
the right is proportional to ip, which in Eq. (|J) leads to 
dissipative effects on ip from its interaction with the x~ 
fields. Analogous expressions also apply to the fermionic 
averages (for Fermi-Dirac statistics) in Eq. (||). 

Based on a systematic perturbative approach, as we 
have shown in previous work we can write Eq. (||), 
using the expressions for the associated averages for the 
Xij, ip and rj fields and Eq. ([?]), as 



V + 3H<p + V^( l p,T) + r ] (v)ip = 



(4) 



where V^ s (ip,T) = dV c ff(tp,T)/d<p, is the field derivative 
of the 1-loop finite temperature effective potential, which 
can be computed by the standard methods (in fact this 
is the resummed effective potential, with masses written 
in terms of the finite temperature effective masses) and 
rj(ip) = r] B (ip) + rp (ip) is a field dependent dissipation, 
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In what follows, we will be interested in the regime where 
r}(ip) > ZH in Eq. (|). As discussed in @,|,@, Eq. (|) 
is valid provided the adiabatic regime for ip is satisfied, 
i.e., the dynamic time-scale for ip must be much larger 
than the typical collision time-scale (~ T -1 ), or \ip/<p\ ^> 
r _1 , where T is the smallest of the thermal decay widths 
r^, r x „, r^, y , as it will set the largest time-scale for 
collisions for the system in interaction with the bath. 

Note that the damping coefficient depends on ip 2 or 
(ip - Mi) 2 . To obtain enough inflation, the amplitude of 
the inflaton should decrease very slowly, which requires 
efficient overdamping. This overdamping is guaranteed 
by having a succession of fields always thermalized, so 
that the population of decay products is not depleted by 
the expansion. This condition is satisfied through our 
choice of Lagrangian for the DM-model, which naturally 
guarantees that a population of decay products will be 
generating efficient damping as ip. slowly rolls down. 

We construct a warm inflation scenario based on the 
following DM-model. For definiteness, we will refer all 
dimensional quantities to Tbi, the temperature of the 
universe at the beginning of warm inflation. The crucial 
property of the DM model of Eq. |l]) is that for a given 
temperature T, only the fields with masses g 2 (ip — Mi) 2 ~ 
T 2 will contribute to the thermal viscosity; the effect of 
heavier fields can be neglected. Thus, as the inflaton rolls 
down its potential, we only must consider the subset of 
fields for values of i which satisfy the above inequality. 
Note that with this choice of model, if ip ^S> Tbi, which 
is needed for efficient inflation, it is still possible to have 
many light Xij (^ij ) fields if ip ~ M, . This will guarantee 
an efficient viscosity term in the equation for ip. 

For convenience, we write the mass parameters as 
Mi = (i + i mia )fh^T BI /g, i = l,...,N M 



a dimensionless constant of 0(1) and M,- 



with m™^ x 

i-mi-aTBl/g 

the lowest mass level coupled to ip. For such a model, 
at T ~ Tbi and Mj+i > ip > Mi, there will be a range 
of masses when 2.5N X + 1 Xij> ?Piji Xi+iJ> V'j+ij'i and r\ 
fields are thermally excited. All other x^^T/v^-fields, 
for i' ^ 1 and j = 1, . . . N x (N x /4), have masses 

™2 



g 2 (v-M,Y > (7 



L x4> 



''Tbi) 2 ~ T and are thus 
10T 



cold (actually thermal excitation begins once m x1 p 
and nontrivial dissipation once m x ^ w (2 — 3)T). 

For the rj — x ~ "P system participating in the dissi- 
pative dynamics of ip, in each interval Mi + i > ip > Mi, 
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the radiation energy density is (2N X + 7N^ +1)tt 2 T 4 /30. 
In addition to these fields, in general there can be a 
set of "heat bath" fields that increase the particle de- 
grees in the radiation system, but do not otherwise con- 
tribute to dissipation. For later use, we addopt the fol- 
lowing notation for the number of these heat bath fields, 
2(15iV x /8+l) 1+Q -(15iV x /4+l), so that, in total, the radi- 
ation system has = 2(15iV x /8 +l) 1+a particle degrees 
for any interval where tp is. a > is a free parameter. 

To simplify the calculation to follow, we consider the 
region where p T and the Xp 3 /6 term dominates the 
equation of motion. From V^ s (p, T) the leading x, ip ne ld 
contribution for M i+1 > ip > M t is (N x T BI /8)g 2 [(<p - 
Mi) + (ip - M;+i)] < 0(gN x T 3 ), so that the constraint 
requires A^ 3 ^> 0(gN x T 3 ). In the perturbative regime 
that we examine, A <C 1, so that m 2 (<yS,T) w Xip 2 /2 = 
fhfyifaTjT^j, since p > T. The xW" masses will 
range from / 2 T 2 /12 (g 2 T 2 /6) < m 2 (m 2 ) < / 2 T 2 /12 

(g 2 T 2 /6) + g 2 (ip - M z ) 2 ~ 0(T BI ). As a simplification, 
the Xj - _masses "vvill be estimated at their largest value 
(m^) 2 = (fh^T BI ) 2 » g 2 (cp — Mi) 2 \ max » 0{T 2 BI ). 
We can then express the condition that the Xp> 3 /Q term 
dominates the equation of motion Eq. (||), in terms of 



"tc<I>/<p — 



3gN x m™*T 2 BI 



4\p 3 BI 



< 1 



(7) 



where, R 



xi>/<p 



_ 3g 2 N x T% I [( v -M i ) + ( v -M i+1 )] I1 



To impose the most stringent constraint from this, it will 
be taken at the maximum value wi™^ x from Eq. (Q). In 
fact, considerable increase in dissipation, thus improve- 
ment in the results to follow, occur by accounting for 
corrections when the x's (^>'s) are in the smaller mass re- 
gion. In the regime outlined above, the effective equation 
of motion for ip, Eq. (0), in the interval Mj+i > ip > Mi 
and in the overdamped limit is 



Vi,i+l{<P)<P ' 



A 3 



(8) 



where = rjf [(ip - Mi) 2 + (ip - M l+1 ) 2 ] + rifT 2 

where rjf and ryf (from and are given by: r)f ~ 
384N x g 4 /[irT(f 2 + 8 3 4 )] \n(2T/m x ) and r?f ~ llN^/T, 
respectively. 

As p moves through the interval Mj + i > ip > Mi, 



TH,i+i{<p) = K(l + r FB )r,f(m™^T BI ) 2 /g 2 , 



(9) 



with r FB = ViT BI g 2 /[r?f n(m^T B i) 2 ) « 0.2, where the 



estimate on the right is from ]9|,|14| for the high temper- 
ature limit and 0.5 < k < 1. Since we are examining the 
region i/5 > T, we can ignore the weak 99 dependence in 
»7t,i+i. In this case, the solution to Eq. (ph is 



<p(t) = PBI 



Wa 



-1/2 



(10) 



with Wi = 3H BlK r,f(m^ T BI ) 2 (l + r FB )/(g 2 \ip 2 BI ) and 
H B j = ^j2Tt\(j) BI /(9m 2 ). Time has been rescaled as t — 

t/H B i with the origin chosen such that ip(r = 0) = ipBii 
which implies tq = 0. 

The resulting warm inflation cosmology for such a field 
trajectory in any interval Mj+i > ip > Mi is equivalent to 
the n — 4 model in |}| . It yields a power-law warm infla- 
tion which never terminates and for which in the steady 
state regime p r (t) / p v (t) = const. In our model, warm 
inflation terminates into a radiation dominated regime 
once ip < Mi = Q, since below that point the dissipative 
term becomes negligible, in which case ip coasts down 
the potential. The essential point of interest here is to 
show that once ip reaches Mo, sufficient e-folds N e have 
occurred while the universe has nonncgligiblc radiation. 

For simplicity, the steady state cosmology in flat spa- 
tial geometry will be examined, which implies from jj), 
for Wi » L p r (0)/Pv(0) = p r (r)/p v (r) = 1/(2W 4 ). 
[From Eq. (|l0|), we can show that this is the necessary 
condition to satisfy the adiabatic condition, \<p/<p\ 3> 
and r^( x ) > H .] In terms of the parameters of the model 
this can also be written as g*w 2 T BI /30 = \p A BI / (ASW4). 
Initial conditions that are more realistic, such as p r (0) > 
p„(0) rapidly evolve into the steady state behavior. In 
this steady state regime, the scale factor solution is B 



R(r) = ( T /W 4 + 1) 



(11) 



with initialization R(0) = 1. 

Our goal is to compute Wi from the microscopic pa- 
rameters of the model and consistent with the many 
constraints given above and in section V of ||. The 
power-law expansion behavior of the scale factor, Eq. 
(|i"T|), is such that the major factor of growth happens for 
r/Wi < 10. N c = Wi e-folds occur at time r/Wi = e-1. 
We will restrict our calculation within this time interval. 
At the end of this interval <p(r), Eq. ([l0|), and T have not 
changed significantly, ip B i — > Pbi/Vg, T bi — > T B i/y/e. 
This simplifies the constraint equations, since they can 
be analyzed at the initial values ipsi, Tbi and approxi- 
mately will be valid throughout this interval. 

The constraint equations for computing Wi are 
as follows. To satisfy the thermalization conditions 
r x ,r0,r</, > H, we will set H B i = min(r x ,r0,r^). 
More explicitly, for the warm inflation solutions studied 
below, r x (which is the smallest of the Ts) is larger than 

Hbi for / ~ 0.8. This condition may be under restrictive 
to obtain good thermalization, but it provides the maxi- 
mum parameter region that may be useful. It should be 
noted that since H ~ ip 2 and r ~ T, as warm inflation 
proceeds, thermalization improves. The thermalization 
condition automatically implies that the adiabatic con- 
dition is comfortably satisfied. 

All the constraints are expressed in the following four 
relations 
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R 



Wa 



15g 4 Y 2 



-(l+a) 



2567r 2 K[min(r ,r x )]r;B(l + r FB ) 
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„max 

x' n xi> 



5127r 2 m 2 i? 2 



R X iP/ v <Pbi/Tbi = 5gN : 



l + Q 



•"£27 (8r 



and 



Ai? 



X</>/V 



= 128?7iS 



(3 3 7V x m-; x )- 



(12) 



(13) 



(14) 



(15) 



Eq. rtl2) is obtained from the relation for W4 below 
Eq. (pp[), where T)(ip) in Eq. (pi), has been expressed 
in terms of R x< j,, Eq. (]?]). The procedure is to input 
/, my, fh^) B j , Y x on the right-hand-side of the above 
four equations, then obtain the left-hand-side of Eqs. 
( |l2| ) and ( |l4| ) from which the remaining two equations 
follow, up to a choice for R x ^/ V < 1. Note that the only 
cases requiring additional heat bath fields (a > 0) are in 
parametric regimes when the right-hand-side of Eq. (|l2] ) 
is greater than 1, since we require R X tp/ v < 1. 

In II analytic expressions were obtained for T(q) and 
?/i in the simplified limit |q| — 0, m x ,m^ ~ 0(m$) 
(which here we call level 1) as well as in terms of the 
exact 2-loop expressions, which must be computed nu- 
merically (level 2). The results in Fig. 1 present both 
levels of approximation. We used g = f = 0.9, k = 0.5, 

' 1 (a ~ 0). The 
; 0.9, for level 1 



N x = 12, = N x /A = 3 and R x ^ /lp 
solid and dashed curves are for to 



and level 2, respectively, and the dotted curve is level 2 
for m™^ x = 2.5. For both curves drawn in Fig. 1, going 
from fhfy BI = 0.002 to 0.05, the initial field displacement 
Vbi/Tbi ranges as 10 6 -10 3 , N M ranges as 10 5 -10 2 and 
A ranges as ~ 10 -17 — 10 -9 . In the regime we are con- 
sidering, i m i„ ~ ipBi I (gTsi) and for the above results, 
imin ranges from 10 2 — 10 6 . To obtain N c ~ 60 e- folds 
of warm inflation for all three cases, fBi/Tsi ~ 3000, 
Nm = 1000, A s=s 10~ 9 . An absolute scale is fixed by set- 
ting to p = 10 19 GeV from which for N c ~ 60 e-folds, we 
findT B/ - (10 13 -10 14 )GeV and H B i - (10 9 -10 10 )GeV. 
[The temperature at the onset of warm inflation, p v = p r , 
is (2W4) 1 / 4 ~ 3.3 times bigger than Tbi, and rapidly de- 
creases to Tbi during the transient period.] 

The DM-model studied here was motivated by the re- 
quirements of warm inflation, dissipative dynamics and 
perturbative renormalizability. We can justify our choice 
of Lagrangian by noting that string-inspired models can 
display an inflaton coupled to mass modes of a string, 
as explained in Ref. jlfl]. Within this context, the large 
number of fields necessary to realize sufficient inflation 
is a natural consequence of the modifications of short- 
distance physics required by string theories. 



In summary, first principles quantum field interactions 
can realize a warm inflation regime with sufficient dura- 
tion to solve the horizon/flatness problems. The inter- 
play between inflationary expansion and radiation pro- 
duction has been a persistent problem since the ear- 
liest history of inflationary cosmology. Thus, despite 
the many questions opened by our model, its underly- 
ing mechanism is a unique and simple resolution to the 
problem. Further study of the inflaton fc-modes is neces- 
sary to address the density fluctuation problem. 
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FIGURE CAPTIONS 

Figure 1: Number of e-folds of warm inflation N e 



versus m^, BI 



for 



0.9 level 1 (solid), 0.9 level 2 



(dashed), and 2.5 level 2 (dotted) with g = f = 0.9, 
K = 0.5, N x = 12, and = N x /A = 3. 
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